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The expansion of a harmonic potential that holds a quantum particle may be realized without any 
final particle excitation but much faster than adiabatically via “shortcuts to adiabaticity” (STA). 

While ideally the process time can be reduced to zero, practical limitations and constraints impose 
minimal finite times for the externally controlled time-dependent frequency protocols. We examine 
the role of different time-averaged energies (total, kinetic, potential, non-adiabatic) and of the 
instantaneous power in characterizing or selecting different protocols. Specifically, we prove a virial 
theorem for STA processes, set minimal energies for specific times or viceversa, and discuss their 
realizability by means of Dirac impulses or otherwise. 
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I. INTRODUCTION 

The need to shorten the times of adiabatic process in 
different quantum or classical systems has motivated a 
surge of activity to design “Shortcuts to adiabaticity” 
(STA), namely, fast protocols for the external parameters 
controlling the system that are generically non-adiabatic 
but provide the same final populations than the adia¬ 
batic dynamics, see [l| and references therein. The time- 
dependent harmonic oscillator is a paradigmatic study 
case, as it represents generically the conditions near equi¬ 
librium and has allowed for experimental implementa¬ 
tions Hi, and applications in many different classical 
or quantum systems 3-H3- For anharmonic corrections 
see as well [3l], [Hj. Since a basic aim of the shortcuts 
is to shorten the process time tf, some “price” in the 
form of high transient energies may be expected to be 
paid. Several articles have studied the energies and pro¬ 
tocol times involved: their characterization; their mu¬ 
tual relation, which is not necessarily of the simple form 
of a time-energy uncertainty principle; and also their 
optimization under different constraints and conditions 
dm nn. In this paper we build on these results and 
continue the analysis of the energies and times in STA 
processes. We do not aim at being comprehensive but 
at examining points that either had not been discussed 
or needed some clarification. Our motivation here is of 
fundamental nature, but the results will be relevant to 
quantify the third principle of thermodynamics i a a, 
as well as the speed, efficiency limits, and costs of quan¬ 
tum engines and refrigerators 3, [3 13 [3 . [32j. 

After a quick review of invariant-based engineering in 
the Introduction, we shall first put forward in Sec. El 
a virial theorem for STA processes which relates time- 
averaged total, kinetic and potential energies, E, K , V. 
We then discuss an idealized protocol that minimizes the 
time-averaged energy for a given time tf, or viceversa, 
by means of Dirac-delta impulses of the spring constant. 
It is shown that these impulses contribute to the total 
averaged energy and we specify how much. After a com¬ 
parison of several protocols in a ( E,tf ) plane we study 


the non-adiabatic energy and its time-average in Sec. IIIII 
Its lower bound for given tf is found and its implementa¬ 
tion is discussed. We also compare different protocols for 
this quantity. Finally, the instantaneous power in STA 
processes is examined in Sec. hyi and we end up with a 
discussion. 


A. Invariant-based inverse engineering 

Here we consider a particle of mass m in a harmonic 
trap, with time-dependent Hamiltonian 

H {t)= 7^ + \muitfq 2 , (1) 

K V 

where ui(t) is the time-dependent angular frequency of 
the trap (from now on “frequency” will be understood 
as an angular frequency), and q and p are position and 
momentum operators, with q being just a c-number in 
coordinate representation. Lewis and Riesenfeld found 
the quadratic invariant [3] I(t) = 7r 2 /2m + mu>Qq 2 /b 2 , 
where u>q = w(0) is the initial frequency, ir = bp — mbq 
is a momentum conjugate to q/b, the dots represent 
derivative with respect to time, and b = b(t) is a time- 
dependent, dimensionless scaling function proportional 
to the width of the elementary dynamical modes. Since 
an invariant I(t) must satisfy —i\I,H\/h + dI/dt = 0, so 
that (d'(t)|/(t)|’I , (i)) remains constant for any |\H(t)) that 
evolves with H(t), b must obey the Ermakov equation 

b + co 2 (t)b=d ( 2 ) 

Any state evolving with H (t) can be expanded as a super¬ 
position of “dynamical modes” l'k(t)) = ]T) n a n \ip n (t)), 
where n = 0,1,2,.., the a n are time-independent am¬ 
plitudes, \ip n (t)) = e l “"(*)| u n (t)), a n (t) = —(n + 
\)wq f 0 dt'/b 2 (t') are Lewis-Riesenfeld phase factors, and 
the | u n {t)) are eigenvectors of the invariant, I(t)\u n (t)} = 
X n \u n {t)), where X n = (n+l/2)Huio are time-independent 
eigenvalues. 
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The dynamical modes have the form 


(glV’nO)) = 


x 


mojQ \ i g-d^+i/ 2 ) fo(ui 0 /b 2 )dt' 

nh ) (2 n n!&) 1 / 2 



(3) 


where T~L n is a Hermite polynomial. If the commuta¬ 
tion relations [66(0),/(0)] = 0 and [66(6/), 7(6/)] = 0 are 
satisfied the initial eigenstates of H are mapped dynami¬ 
cally into final eigenstates, avoiding final particle excita¬ 
tion These commutation relations indeed hold if the 
boundary conditions 


6 ( 0 ) = 1 , 6 ( 0 ) = 0 , 

6 (tf) = 7, b(tf) = 0, 


are satisfied, where 7 = ^uio/ojf , and w/ is the final 
frequency. In addition, 

6(0) = 6 ( 6 /) = 0 (5) 

may be imposed to make w( 6 ) continuous at 6 = 0 and 
6 = tf. To inverse engineer w( 6 ), 6 ( 6 ) is designed first, 
interpolating with some convenient function between the 
boundary conditions at 0 and 6 /, and w( 6 ) is deduced 
from the Ermakov equation © Q. A simple protocol 
based on a quintic polynomial for 6 ( 6 ), where the coeffi¬ 
cients are found from Eqs. © and ©i is 0 

6 ( 6 ) = 6(7 — l)s 5 — 15(7 — l)s 4 + 10(7 — l)s 3 + 1 , ( 6 ) 

where s = t/tf. For very short times, imaginary frequen¬ 
cies appear for some intermediate times, corresponding 
to a negative, concave-down potential. 

Another simple solution is the “bang-bang” (stepwise 
constant) form of the control function ui(t). It results in 
particular from applying optimal control theory (OCT) 
to minimize the time for given constraints on the fre¬ 
quency 0 , a in Hi- Some relevant expressions are pro¬ 
vided in the Appendix. As well, “bang-singular-bang” 
solutions result from minimizing the time-averaged en¬ 
ergy for the same constraints [l2|. 


II. TOTAL, KINETIC AND POTENTIAL 
ENERGIES 


For the n-th dynamical mode, the instantaneous en¬ 
ergy, E n (t) = ( ip n \H(t)\ip n ), is given by 


£n(6) = 


(2 n + 1)6, 

4w 0 


6 2 + W 2 (6)6 2 + ^ 
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(7) 


We may divide the energy into kinetic and potential parts 
using the separation H = K + V in Eq. ©. Then 
E n (t) = K n (t) + V n (t), where 


K n {t) 

V n {t) 


(2n + 1)6 /• 2 wg 

4co 0 \ 6 2 


(2n + 1)6 

4wq 


6 2 w 2 (6). 


( 8 ) 

(9) 


The mean value theorem implies for the interval (0,6/), 
assuming that the boundary conditions © are satisfied 


bmax > \bmax\ > 2^1. (10) 

6/ t f 

The first inequality sets a lower bound for the max¬ 
imum of the instantaneous kinetic energy, K n ^ max > 

f / ) . In general we cannot set bounds for mini¬ 
mal or maximal values of V n unless physically motivated 
bounds apply for both 6 and w. 


A. Virial theorem for STA processes 

A virial-thcorem relation will be proved for STA pro¬ 
cesses that implies, remarkably, that the time-averaged 
energy is equipartitioned into kinetic and potential con¬ 
tributions, namely, K n = V n = E n / 2, where the over line 
denotes a time average over [0,6/], i.e., for a generic A n (t) 
function, 


A n = [ A n (t) dt. (11) 

tf Jo 

The virial relation for STA processes is not at all an ob¬ 
vious result, as the ordinary quantum-mechanical virial 
theorem applies to localized (square integrable) station¬ 
ary waves and a time-independent Hamiltonian |3fil - l38j . 
Generalizations were proposed for time-dependent states 
and a time-independent Hamiltonian, performing the av¬ 
erage over an infinite time or over a period if the motion 
is periodic [tI [© . Finally, in Ref. [42] a virial theorem 
was found for a charged particle in a time-dependent elec¬ 
tromagnetic field. 

Suppose first that the trap expansion is performed adi- 
abatically over a very long time 6/ so that the state re¬ 
mains in level n from 6 = 0 to tf. The ordinary virial 
theorem, i.e., the one formulated for stationary states, 
can then be applied. For the harmonic oscillator in level 
n , it states that V n = K n = E n /2. Along the adiabatic 
process the values of the energies will change slowly in 
such a way that the relation is preserved. We may thus 
take a time average and get V n = K n = E n / 2. For an ar¬ 
bitrary non-adiabatic process, however, this relation will 
not be true in general. 

According to Eq. © and Eq. (fTTT) . the time-averaged 
energy is 


1 f tf (2n + 1)6 


tf 


4 w 0 


6 2 + w 2 (6)6 2 + -A dt, (12) 
b z 


or, using Ermakov’s equation ©, 

E n = — 

tf 


f tf (2 n + 1)6 

2^ + 6 2 - 66 

Jo 4w 0 

[ 6 2 


dt. (13) 
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By using the boundary conditions Q and partial inte¬ 
gration, we have f Q f (-bb)dt = f 0 f b 2 dt. Then the time- 
averaged energy can be rewritten as, see Eq. flHJ, 


E n — E n o — 


1 r t} (2n + l)h 
tf Jo 2 wq 


b 2 


b 2 


= 2K n . (14) 


(We have introduced a specific notation, E ni2 for the in¬ 
tegral in Eq. m to emphasize that E n ^2 is only the 
averaged energy when the boundary conditions are sat¬ 
isfied.) Since E n = K n + V n , we find, remarkably, the 
virial-theorem relation 


K n =V n = E n / 2. (15) 

This implies a number of consequences. For example, 
minimizing any of the time-averaged energies automati¬ 
cally minimizes the others. As well, non-trivial bounds 
may be set since bounds for one of the averaged energies 
works for the others as well. Thus, in spite of the fact 
that V n may be negative during some time interval, the 
time-average V n must be positive, as K n > 0; even more, 
it will be bounded from below by a positive number as 
described in the next subsection. Further applications 
in this work may be found in Sec. IIII Al and in the Ap¬ 
pendix. 


B. The energy contribution from Dirac impulses 

Using the Euler-Lagrange equation [§| to minimize Eq. 
(fT4l) or Optimal Control Theory (OCT) [l^, the follow¬ 
ing quasi-optimal function is found, 


b(t) = (B 2 - u%tf 2 ) (J-^J + 2b(±)+1, (16) 


where B = to^t 2 + y 2 — 1 and the positive root should 
be taken; this only satisfies the boundary conditions 
&( 0 ) = 1 and b(tf) = 7 at the time limits but not the 
conditions on the derivatives. Therefore, as such it only 
provides a bound for the time-averaged energy. Substi¬ 
tuting Eq. © into Eq. ©. this lower bound is 


when the cap duration r went to zero, E n —> E n l. How¬ 
ever the caps did contribute to the integral, actually half 
of the total time-averaged energy came from them, but 
the significance of this fact was not discussed. 


In 


12 ] an alternative to the polynomial caps was put 


forward, namely, Dirac-delta impulses of the control func¬ 
tion u> 2 (t) that switch the derivatives 6 ( 0 + ) and b(tj) to 

b(0~) = b(tfi = 0, 


w 


t < 0 “ 


D 0 S(t), 0 - < t < 0 + 


\t) = 


D f S(t) 


LOf 


l 0 +<t<K 

tj <t <tf 

/+ 


(18) 


tj<t 


Substituting the (5-terms in the Ermakov equation 
integrating b around (in the immediate neighborhood of) 
t = 0 , and tf, and taking into account the boundary 
conditions in (H]), D 0 and Df are found to be 


D 0 = - 


6(0+) 

m 


, _ 

f b{t f ) ■ 


(19) 


(The first impulse always corresponds to a negative 
Delta, while the second may have the two signs 0 .) The 
protocol (1181) with Dirac impulses provides a formally el¬ 
egant proof that the bound can indeed be reached, at 
least in principle, because, as the boundary conditions 
are satisfied, Eq. (fldl) may be used to get Eq. (El as 
the actual time-averaged energy of the process. Note 
that these protocols minimize the time-averaged energy 
for a given time tf, but also minimize the time tf for 
a given time-averaged energy. In the idealized processes 
that realize the bound the instantaneous potential energy 
jumps to (plus or minus) infinity due to the Dirac pulses. 
It might seem that the deltas do not contribute to the 
averaged energy since uj does not appear explicitly in Eq. 
m- However they do. To see why and how much, the 
time-averaged energy given by the original expression, 
Eq. (1121) . where no partial integration has been carried 
out, can be separated into two parts, one from the deltas 
and one from the central time interval, 


E n =A s + E n [0 + ,tj], ( 20 ) 


EnL = 


(2 n + 1)H 
2u>ot 2 


(B 2 - Lu^t 2 ) - 2 uj 0 tf 


( B 2 + B - oj 2 t 2 f \ 

arctanh -— — arctanh 

V w 0 tf 


B 

Wo tf 


where 


(17) 


A 5 



E n (f) dt + 



E n (t) dt, 


E n [0 + ,tj] 



E n (t) dt, 


( 21 ) 

( 22 ) 


which becomes E u l ps ^ 2 ^^ f° r tf *C 1 /^/woW/ and 

7 1. In 3, polynomial “caps” were added around a 

central time-segment defined by Eq. EK to match this 
function with the right boundary conditions. A numeri¬ 
cal example for tf <C 1 /^/wqw/ and 7 1 showed that 


and E n (t) is given in Eq. (0 with the frequencies in Eq. 
m- Making use of the coefficients m, the contribution 
from the Dirac impulses to the time-averaged energy is 


Ag 


(2 n + l)h 
4w 0 tf 


HtfMtf) - H 0 + ) 6 ( 0 ) 


(23) 
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The other term in Eq. (1201) can be rewritten using partial 
integration as 


E n [0^~ 1 t j] — En, 2 T A boundary ; (24) 

where 


A 


boundary 


(2n + 1)H 
4 uj Q tf 


Thus it turns out that 


~ b{0 + )b(0) . 

(25) 


A,5 — A boundary i 

so, according to Eqs. 


E n — E n 2 — E u l. 


(26) 

(27) 


Substituting the quasi-optimal protocol (fT7>l) in Eq. (1^1) , 


A s 


(2 n + 1 )h 
4 coot} 


(B 2 



When tf <C 1/^UJquJJ and 7 > 1, 


A^ 


(2n + l)/i 

4w/^ 


Eul/ 2 , 


(28) 


(29) 


which agrees with the contribution of the caps found nu¬ 
merically in [&] as t — \ 0. The reason for this result did 
however escape the authors of Q. 


C. Comparison of protocols 

Figure Q] depicts tf versus the time-averaged energy for 
different protocols. A log-log representation is chosen to 
show the global behavior and the domains where the pro¬ 
tocols are applicable. As announced the fastest protocols 
for a given E n are the ones that implement the bound 
dm. Other constraints lead to different winners. For 
example, if |w 2 | is limited by some predetermined value 
the fastest solutions are of bang-bang form 0 , 0 . If in 
addition to this bound the time-averaged energy is fixed, 
they are of bang-singular-bang form [121 ] , where the pro¬ 
tocol in Eq. dm applies in a central segment flanked 
by constant-frequency intervals. The bang-bang proto¬ 
col considered in Fig. Q] assumes the two intermediate 
frequency steps iuj\ and wi for durations t\ and £ 2 , re¬ 
spectively; wi decreases for larger times, and a maximal 
time (minimal E n ) exists, as discussed in the Appendix 
and shown in the figure. 


III. NON-ADIABATIC ENERGY 

The “non-adiabatic energy” is defined as the difference 
between the total energy and the energy of a correspond¬ 
ing adiabatic process, 

E na (t) = (H(t)} - E ad (t) 

= E[^(*)-^(°)]^)’ (30) 

3 



log I0 [EJE,m 

FIG. 1: (Color online) Protocol time tf versus time-averaged 
energy E n for a quintic polynomial protocol (dotted brown), 
bang-bang for uo\ = W 2 , which decreases for larger times, see 
Eq. (EH (solid black), and lower bound (dashed red). The 
filled black point is for the bang-bang corresponding to — to 2 < 
u 2 (t ) < to 2 . The bang-singular-bang form (dot-dashed blue) 
is for the same bounds, with the bang-bang control of the 
point as its limiting, minimal-time case. Parameters: to 0 = 
2500 x 2-7T Hz and to/ = 25 x 2n Hz. 


where the £j(t) = ( j + 1/2) hio (t) are the instantaneous 
eigenenergies, the Vj(t) the corresponding populations, 
and E a d = JT Vj(0)ej(t). We must now assume w(f) > 0 
in order to have a meaningful real quantity. The “mini¬ 
mal work principle” [0], establishes that ( H{t)) > E a d(t) 
provided that the initial state is passive (i.e., the initial 
density matrix is diagonal in the energy representation 
and satisfies V n {0) > V n +i(0)) and no level-crossings oc¬ 
cur. We shall restrict this section to the ground state 
|^o) as a simple passive state. Then, from Eq. flTJ), 

E na (t) = ^- 0 ( i}2 + “ W + - ^Mt). (31) 

In the isentropic expansion stroke of an Otto cycle the 
non-adiabatic energy gives the dissipated work, that 
must vanish at tf for a STA process. Its time average 
has been proposed to quantify the cost to implement a 
quantum engine based on the harmonic oscillator [27]. 

The time-averaged non-adiabatic energy is 

_ = ^o [** /E.EM . i_mEu 
4 tf Jo wg b 2 w 0 J 

(32) 

Using Ermakov’s equation, partial integration and the 
conditions a, we get a simpler form, 


E na = E% a 


fvjJo j' tf ( 

Jo y w o 


1 



(33) 


where, again, we have introduced a special notation for 
the last integral, E J", to emphasize that it becomes the 
time average of the non-adiabatic energy provided the 
boundary conditions a are satisfied. 
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FIG. 2: (Color online) A classical particle (blue dot) 
moves along the bottom of an expanding potential U(t) = 
[cj(t) 2 b 2 + Wq/ 6 2 ] /2 (red line), see text. The two parts of 
17(f) are shown, uj(t) 2 b 2 /2 (blue dashed line), and wo/( 2 b 2 ) 
(black dashed line). 


A. Lower bound for time-averaged non-adiabatic 
energy 

A useful analogy exists between the Ermakov equa¬ 
tion and the dynamical equation of a fictitious (classical!) 
particle with mass m = 1 and dimensionless “position” b 
moving in a potential of the form (44j 

U(t) = i [u 2 (t)b 2 +Wg/ 6 2 ] (34) 

with kinetic energy 6 2 /2. Newton’s equation takes indeed 
the form of the Ermakov equation, 

6 = Wg/ 6 3 — w 2 (t) 6 . (35) 

This analogy provides useful intuition and will help us 
to find the lower bound for E na . The corresponding (di¬ 
mensionless) energy is 

■* = i(i a + ^(<)!’ 2 + )f). (36) 

Up to a constant factor, it has the same form than the to¬ 
tal energy 0, Eq = Beware that the kinetic and 

potential energies of the classical particle do not corre¬ 
spond in general to the quantum counterparts in Eqs. 
<18191) . Specifically, the term Wq/(2 b 2 ) is “kinetic” in the 
quantum scenario, and “potential” in the classical anal¬ 
ogy. 

Let us now calculate the excitation energy of the par¬ 
ticle E ex measured from the potential minimum. To find 
the minimum we calculate dU{t)/db = 0, which leads to 

“<*> = 5 % ■ < 37 > 

and U m i n = loqui, so that 

E ex = 1^ + 1 (“> 2b2 + ^ - 2w 0 w) • (38) 


This has exactly the same form than the non-adiabatic 
energy (1311) up to the scaling factor, E na = ^ E ex . To 
complete the analogy, we consider that oj{t) changes from 
ui o to ujf and trajectories from the initial to the final 
points of the potential minimum, 5(0) = 1 to 6 (t/) = 7 . 

Note that both the first term (kinetic energy) and the 
second one (potential energy measured from the mini¬ 
mum) are positive. The potential part of Eq. (1381) can 
in fact be made zero if the particle moves all the time at 
the bottom of the potential, without ever being affected 
by a force, then Eq. (l37l) should hold and, substituted in 
Eq. (l35l) this gives the equation 

6 = 0, (39) 

with solution 

6=1 + (40) 

*f 

that satisfies the boundary conditions 6 ( 0 ) = 1 and 
b(tf) = 7 . This linear 6 also minimizes the time-averaged 
kinetic energy for these boundary conditions, as it can 
be seen from the Euler-Lagrange equation. However the 
particle trajectory described by Eq. (l40l) is not at rest at 
t = 0 and tf, so for the more restricted family of trajec¬ 
tories satisfying Eq. 0 it only provides a lower bound 
for E ex . Of course the same (scaled) lower bound is valid 
for the analogous quantum system, namely 


^na _ 


h 

4 Wo 



(41) 


When 7 > 1, if£ a « ^ , which is half the lower 

f f 

bound of the time-averaged energy for the ground dy- 
namical mode if, in addition, tf -C 1 //. Under 
these conditions the quantum kinetic energy is dominated 
by the first term in Eq. ©, which corresponds to the 
classical kinetic energy, and the virial theorem implies 
K~ 0 = E~ 0 / 2. 


B. Comparison of protocols 

The linear trajectory of 6 ( 6 ) in Eq. (UU1) does not sat¬ 
isfy all the boundary conditions 0 and thus the protocol 
ui(t) = coo/b 2 based on it leads to excitation. Dirac im¬ 
pulses are not an option now, as the one at t = 0 would 
have to be negative but uj 2 < 0 is not allowed. We may 
instead complement the protocol with caps of durations 
tl and rg (for “launching” and “stopping” respectively) 
that connect the linear function with the proper bound¬ 
ary conditions, similarly to Q. A hybrid protocol defined 
in this manner is 

f En=0 friS n , 0 <t<T L , 

= < (7 “ i ) 3 + 1 , t l <t <t f - t s , 
lEn=0 2nS n , tf -T S <t<tf, 


(42) 
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FIG. 3: (Color online) Dependence oft/ on the time-averaged 
nonadiabatic energy for the hybrid protocol in Eq. (1421) (dot¬ 
ted brown line), quintic polynomial protocol (solid blue line), 
bound in Eq. (1411) (dashed black line), and bang-bang proto¬ 
col with different bounds (dot-dashed red line). The param¬ 
eters are wo = 2500 x 2n Hz and w/ = 25 x 27r Hz. 


where s = t/tf, the coefficients f n are found from the 
equations that match 6 and & at 0 and t£ , and the g n from 
the matching at tf — ts and tf. It can be proved that 
in the cubic interpolation b > 0. The cap durations tl,s 
can be adjusted by a subroutine that minimizes E na with 
the constraint w > 0. Due to the constraint the cap times 
cannot be made zero, so, unlike Sec. Ill Bl the bound is 
not reached. The caps can only be constructed in this 
way until a minimal tf for which imaginary frequencies 
appear. 

We also consider bang-bang protocols that provide a 
meaningful non-adiabatic energy with wi = 0 , and W 2 = 
ujq/3, /3 > 0, see the Appendix, with 


t i 


t2 


'(7 2 -l)(7 2 /? 2 -l) 

2/32 




1 

-- arcsm 

LUqP 



(43) 

(44) 


Fig. [3] shows the scaling of tf and the time-averaged 
non-adiabatic energy. Note that there is some comple¬ 
mentarity among the different protocols, which are ap¬ 
plicable in different domains. Whereas the bang-bang 
protocol only satisfies the boundary conditions in a small¬ 
time, high-energy regime, the quintic polynomial dH]) or 
hybrid protocols apply rather in a small-energy, low- 
energy scenario. 


IV. POWER 

We finally consider the instantaneous power during the 
STA process, 


P n {t) = 


dE n 

dt. 


k n tb 2 , 


(45) 
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FIG. 4: (Color online) Relative power versus s = t/tf for the 
polynomial ((Gj) (dashed red line) and for the polynomial (1501) 
(solid blue) with C 3 = 78.5088, C 4 = —459.7638. Parameters: 
wo = 2500 x 2-7T Hz, LOf = 25 x 27r Hz, and tf — 8 ms. 


where we have used Eq. (0 and defined k n = 4 " o ’ 

and /C(f) = w 2 . This is a remarkably simple expression 
for P n . (Defining the time-independent operator v by 
V = K.(t)v, and using Eq. ([9|) , it may also be written 
as JC(v) which is the form given in [45J.) For any STA 
process its integral is 

dt = E n (tf) - E n { 0) = ^ h(ujf - w 0 ). 

(46) 

As an example of the possible use of power for shortcut 
design, suppose that we are interested in minimizing the 
power peak (maximum) of |P„|, distributing the power 
homogeneously from 0 to tf to facilitate the power ex¬ 
traction. We then set 

P n = k n tb 2 = C n , (47) 

where C n is a constant adjusted to satisfy Eq. 60, 
namely 



C n = (n + l/2)h(w f -u 0 )/tf, (48) 

note that C n /k n does not depend on n, so that the rel¬ 
ative power P re i = P n /C n does not depend on n. Ac¬ 
cording to the mean value theorem, 1 is a lower bound 
for the maximum of P re i in an arbitrary STA expansion 
process. /C and b must satisfy the Ermakov equation ([2]) 
so Eq. S3 becomes 

2 

b b — bb + 4^-6 = 2(wo — w/)^. (49) 

This is a third order equation. To avoid a sudden jump 
in w(t) at time t = 0, and thus a Dirac delta in the power, 
we must impose 6 ( 0 ) = 0 so in fact, with 6 ( 0 ) = 1 and 
6 ( 0 ) = 0 we already fix the three possible constants at 
time t = 0. The numerical solution may satisfy b(tf) = 7 
for a specific tf but the conditions b(t /) = b(t /) = 0 fail in 
general. We may thus resort to more modest objectives. 
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As an example of many possible strategies, in Fig. [4] 
we plot P re i for the fifth order polynomial (0 and for a 
polynomial that satisfies the boundary conditions (j4| and 
([5]) with two free parameters that are chosen to minimize 
numerically the maximum of P re u 

b(t ) = 1 +C 3 S 3 +C 4 S 4 —(21+6c3+3c4— 217 )s 5 
+ (35+8c3+3c4—357)s 6 —(15+3C3+C4—157)s 7 , (50) 

where s = t/tf. 


V. DISCUSSION 

Different transient energies in shortcuts to adiabatic 
harmonic expansions for a single particle have been stud¬ 
ied, as well as the relations with the process time. This 
is important to establish operational limits, which are 
different from the naive application of a time-energy un¬ 
certainty relation. General bounds are provided, and 
protocols that realize or approach them have been dis¬ 
cussed. The specific experimental conditions, i.e. the 
actual realization of the system and the potential by dif¬ 
ferent interactions (optical, magnetic, electrostatic), will 
determine their realizability, since optimal protocols of¬ 
ten imply Dirac deltas (e.g. in the potential, or in the 
power). 

The virial theorem for time averaged energies emerges 
as an important relation for STA processes. We have 
given examples where it is instrumental in providing 
bounds and interpreting the results. Further relations 
among various energies and their applications will be dis¬ 
cussed elsewhere, extending the theorem to other opera¬ 
tions and potentials. 

The power in STA has been also brought to the fore. 
It provides a further criterion to choose among different 
shortcuts. Even though its time integral is invariant for 
all STA, the maximal values could be minimized. 
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Appendix A: Bang-bang protocols 

A “bang-bang” protocol is composed by time segments 
with constant frequency so that the energy rate transfer 
(power) from or to the system consists of Dirac deltas at 


the switching times. For applications and discussion of 
their relation to Optimal Control Theory see s a mi, 
HQ. Typically the frequency is chosen at the extreme 
values allowed to achieve optimal results. 

We consider here processes with two intermediate steps 
of the form (more switching times may be required to find 
the true minimal time in general El) 


w(t) 


wo, t — 0 

iui i, 0 < t < ti 

W 2 , tl < t < tl + t 2 ' 

W/, t = tf = tl + t 2 


(Al) 


where wi > 0 and w 2 > 0 are real constants, and t\ and 
t 2 are the durations of each part. For 0 < t < t\, the 
solution of the Ermakov equation satisfying 6(0) = 1, 
and 6 ( 0 ) = 0 is 


b{t) = 4 1 + — sinh 2 (wit), (A2) 

V w i 

and for t\ < t < tf, the one that satisfies b(tf) = 7 , 
b(t f ) = 0 is 


b(t) 



w o - 7 4 w| 
7 2 wf 


sin 2 [w 2 (t/ - t)]. 


(A3) 


ti = —arcsinh, 

Wi 1 


t 2 = — arcsin , 
w 2 


M(7 2 -l)(7 2 wl—wg) 

/ 7 2 (w 2 +w2)(w 0 2 +w2) ’ 

^( 7 2 -l)( 7 M+^) 

(w 2 +w 2 )(7 4 w 2 -w 2 ) ' 


(A4) 


There is an upper bound for tf. According to Eq. (IA4I) . 
w 2 > y/uoWf must be satisfied to make t\ > 0. With 
this condition, t 2 will be positive. At w 2 = lJqlJJ , t\ = 
0 , which gives the upper upper bound for tf and the 
corresponding minimal time-averaged energy, 


imax 

v 


E™ n 


7T 


2-y/WoW/ 

(2n + l)H^l. 


(A5) 


To calculate E n , we use the fact that the total energy re¬ 
mains constant during the constant-frequency intervals. 
It can be calculated with Eqs. ([HI) and m, at t = 0 + for 
the first segment and at t = tj for the second segment, 


E n 


(n + 1 / 2)6 

2 


Eo-w^/w 0 , 0 < * < tl, . 

(ej+<4)/cj f , ti<t<t f . ( Ab > 


If u)\ = w 2 wo with 7 > 1, tf approaches zero, and 
the corresponding time-averaged energy is 


Er, 


(2 n + 1)7 tH 
16 w/t 2 


In [ 2 7 ], 


(A7) 


















This is of course larger than the asymptotic value of lower 
bound for the time-averaged energy, E u l ~ . 

tCUftf 

If we set uj\ = 0, 0 J 2 wo, and 7 1 then 

1 

tr S3 t\ S3 , 

7 (A8) 

E n ss (n + l/2)fiw 0 . 

or, making use of the relation between and wo in 
Eq. (IA8I) we find again the bound value, namely E n ss 


^ 2 ■ We ma y interpret this result in view of the virial 

theorem. Under these conditions most of the protocol 
time is just a free expansion ended up by a short seg¬ 
ment with a strong potential. The frequency switch at 
time zero suddenly turns off the initial harmonic poten¬ 
tial so that half the initial energy vanishes, and for most 
of the time the energy (purely kinetic) is E n ( 0)/2. Thus 
K n ss E n ( 0)/2, and according to the virial theorem the 
last segment must be such that V n = E n { 0)/2 as well. 
That explains why E n S 3 E n ( 0). 
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